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We investigate various possible definitions of an effective temperature for a particularly simple 
nonequilibrium stationary system, namely a heated Brownian particle suspended in a fluid. The 
effective temperature based on the fluctuation dissipation ratio depends on the time scale under 
consideration, so that a simple Langevin description of the heated particle is impossible. The short 
and long time limits of this effective temperature are shown to be consistent with the temperatures 
estimated from the kinetic energy and Einstein relation, respectively. The fluctuation theorem 
provides still another definition of the temperature, which is shown to coincide with the short time 
value of the fluctuation dissipation ratio. 
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In the recent years, so called "active colloids", i.e. col- 
loidal particles that exchange with their surroundings in 
a non Brownian manner, have attracted considerable at- 
tention from the statistical physics community [1]. These 
systems are of interest as possible models for simple liv- 
ing organisms, and the description of the corresponding 
nonequilibrium states using the tools of standard statis- 
tical physics raises a number of fundamental questions 
[2, 3]. The most widely studied active colloids are those 
that exchange momentum with the supporting solvent in 
a non stochastic way, resulting into self propulsion. A less 
studied possibility is that the colloid acts as a local heat 
source and is constantly surrounded by a temperature 
gradient. Experimentally [4], such a situation is achieved 
when colloids are selectively heated by an external source 
of radiation which is not absorbed by the solvent. If the 
heat is removed far away from the particle, or, more prac- 
tically, if the particle concentration is small enough that 
the suspending fluid can be considered as a thermostat, 
a simple nonequilibrium steady state is achieved. Each 
colloidal particle is surrounded by a spherically symmet- 
ric halo of hot fluid, and diffuses in an a priori Brownian 
manner. The diffusion constant of such heated Brownian 
particles was experimentally shown to be increased com- 
pared to the one observed at equilibrium [4] , and a semi 
quantitative analysis of this enhancement was presented 
in reference [5] , based on an analysis of the temperature 
dependence of the viscosity. 

In this report, we use simulation to investigate in de- 
tail the statistical physics of the simple non equilibrium 
steady state (NESS) formed by a heated particle sus- 
pended in a fluid. The most natural way of describing 
such a system, in which the particles diffuse isotropically 
in the surrounding fluid, is to make use of a Langevin 
type equation for the center of mass velocity U, involv- 
ing in general a memory kernel and a random force 
R(t): 

= - f ((t- s)U(s) ds + F cxt + R(t). (1) 



In a system at thermal equilibrium at temperature T, the 
correlations in the random force and the friction kernel 
are related by the standard fluctuation dissipation the- 
orem, (R a (t)Rp(t')) = 6 a p({\t - t'\)k B T [6]. Obviously 
such a description is not expected to hold for a heated 
particle, as the system is now out of equilibrium. A gen- 
eralization of Eq. 1, involving a corrected fluctuation 
dissipation relation with an effective temperature T c g re- 
placing the equilibrium one, would however appear as a 
natural hypothesis. In fact, such an approach was shown 
to hold for sheared systems kept at a constant tempera- 
ture by a uniform thermostat [7], or in the frame of the 
particle for a particle driven at constant average speed [8] . 
The interpretation of recent experiments [3] also makes 
implicitly use of such a description in describing the sedi- 
mentation equilibrium of active particles, or in analyzing 
the diffusion constant for hot Brownian motion [5] . 

The use of a Langevin equation with an effective tem- 
perature has several direct consequences. The kinetic 
energy associated with the center of mass, (^MU 2 ), is 
necessarily equal to the effective temperature §fcBT c ff. 
The diffusion coefficient D and the mobility under the 
influence of an external force fi = U x /F x are related by 
an Einstein relation, D//i — k&T e f[ [9], More generally, 
this relation can be seen as the steady state version of 
the proportionality between the time dependent response 
function to an external force, x(t) = SU x (t)/5F x , and 
the velocity autocorrelation in the nonequilibrium steady 
state: 



x(t) 



1 



eff 



-(U x (0)U x (t)). 



(2) 



This relation was explored numerically for self propelled 
particles in reference [2] , and shown to be consistent with 
the observed Einstein like relation. Independently of the 
use of a specific Langevin model, this relation defines 
an effective temperature trough a so called "fluctuation 
dissipation ratio" . The applicability of an effective tem- 
perature description is determined by the dependence of 
this fluctuation dissipation ratio on time. We show in 
the following that the time scale at which the fluctuation 
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FIG. 1. Left- Snapshot of the simulated system for T p = 
3.5e/fcB (To = 0.75e/fcB); Gray levels indicate the kinetic en- 
ergy of atoms. Right- Steady radial temperature and density 
profiles for this system. 

si dissipation ratio of a heated particle is determined indeed 

82 matters, so that a single temperature description, even 

83 in such a seemingly simple system, is problematic. 

84 Finally, the use of a Langevin description with an ef- 

85 fective temperature entails the validity of several "fluc- 

86 tuation relations" [10], which have been the object of 
s? numerous recent experimental and numerical tests, both 
as in equilibrium and nonequilibrium systems. The study of 

89 the fluctuation relation for the heated particle constitutes 

90 the last part of this report. 

91 Our work is based on a direct molecular simulation 

92 (MD) approach of a crystalline nanoparticle diffusing 

93 in a liquid. The simulation were carried out using the 

94 LAMMPS package [11]. Details of the model can be 

95 found in previous works [12, 13], where we used this 

96 system to investigate heat transfer from nanoparticles. 

97 The particle was made of 555 atoms with a fee struc- 

98 ture, tied together using FENE bonds. The liquid was 

100 made of <~ 23000 atoms (Fig. 1). All liquid and solid 

101 atoms interacted via the same Lennard- Jones (LJ) po- 

102 tential v = 4s[(a/r) 12 — (cr/r) 6 ], at the exclusion of solid 

103 atoms directly bonded to each others. In the following, 

104 all results will be given in LJ units, namely a, e/k^ and 

105 t = \J ma 2 /e for length, temperature and time, respec- 

106 tively. The atoms in the solid particle were held at con- 

107 stant temperature T p using a Nose-Hoover thermostat, 
los after subtracting the velocity of the center of mass. In 
109 order to mimic the bulk liquid - far from the particle 
no - acting as a thermal bath, a rescaling thermostat was 
m applied only to liquid atoms lying beyond 15cr from the 
U2 center of the particle (this condition being evaluated each 
n3 time the thermostat was applied), to keep them at con- 
u4 stant temperature To = Q.Toe/k^. This amounts to an 
us assumption that the temperature profile around the par- 
H6 tide follows the latter instantaneously. This is a reason- 
n7 able assumption, as heat diffusion is much faster than 
us mass diffusion in our system: £>hcat ~ la 2 /r [13], while 
u, fl mass G [0.002; 0.02]cr 2 / T (Fig- 2). Finally the whole 
120 system was kept at fixed pressure p = 0.0015e/cr 3 us- 



121 ing a Nose-Hoover barostat. Simulations were run over 

122 typically 10 7 timesteps in order to accumulate enough 

123 statistics. 

124 In previous work, we have shown that nanoparticles 

125 are able to sustain extremely high heat fluxes, via two 

126 mechanisms: Firstly, interfacial thermal resistance at the 

127 nanoscale generates significant temperature jumps at the 

128 interface, i.e. the contact temperature T c of the liquid at 

129 the nanoparticle surface is much lower than the particle 

130 temperature T p (Fig. 1). Secondly, the large curvature- 

131 induced Laplace pressure prevents the formation of a va- 

132 por layer at the interface; At the highest temperatures, 

133 only a stable depleted region is observed (Fig. 1). 

134 Two approaches were used to measure the effective 

135 temperature of the particle. We started by measuring 

136 the kinetic temperature Tk, related to the center of mass 

137 velocity of the particle. Due to the finite ratio between 

138 solid and liquid masses, care has to be taken to measure 

139 the relative velocity between the solid nanoparticle and 

140 the liquid Ui = U S i — Uu (i = x, y, z), with U S i and Uu the 

141 velocities of the solid and liquid centers of mass along the 
.42 i direction. Tk was then given by ^/cbTk = ^ m ctf(U?}, 
143 where m e e = m s m\/ (m s + mi) [m s and m\ being the to- 
rn tal mass of the solid and liquid components] . We checked 

145 that this procedure behaved correctly for all mass ratios, 

146 even when the mass of solid atoms is increased artificially 

147 up to the point where m s = mi. All the velocity measure- 
rs ments presented in the following were done consistently 

149 using this procedure. Tk was evaluated along the 3 de- 

150 grees of freedom of the particle in order to estimate the 

151 uncertainties, which were below 1%. 

152 We also measured the "Einstein" temperature Te, de- 

153 fined as the ratio between the diffusion coefficient D and 

154 the mobility /j, of the particle [9]. The diffusion coeffi- 

155 cient was computed as the plateau value of the integral of 

156 the velocity autocorrelation function (VACF) Cuu{t) = 

157 (Ui(t) Ui(0)) of the nanoparticle: D = lim^oo T>(t), with 
iss V(t) = J Q Cjju(s)ds (Fig. 2a). The plateau is reached 

159 after a correlation time typically around t c ~ 30r. The 

160 mobility /i was computed by applying an external force 

161 F = I0e/a to the particle, and measuring its steady ve- 

162 locity U in the direction of the force: fj, = U/F (linear 

163 response in the applied force was carefully checked). 

164 In Fig. 2.b, we have plotted both measures of the par- 

166 tide's effective temperature as a function of T p . One 

167 can note that all temperature estimates collapse to To at 

168 equilibrium. A striking feature of Fig. 2.b is that the two 

169 approaches to measure the effective temperature of the 

170 particle provide different results. While this is expected 
(7i for active colloids with a ballistic motion at short times 

172 [3] , it is quite surprising in the case of a simple Brownian 

173 particle, and cannot be understood in the framework of 

174 a Langevin description. As discussed before [5] , one can 

175 finally note that neither Tk nor Te identify with the con- 

176 tact temperature T c , as could be naively expected [14] 

177 (Fig. 2.b). 
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To understand the existence of two temperatures in 
the system, we have probed the fluctuation dissipation 
theorem (FDT) for the Brownian system under study. 
Generally speaking, considering a physical observable A, 
the response of a system driven out of equilibrium at time 
t = by the action of a small external field J-(t) is char- 
acterized by the susceptibility XAc{t) = ^fo) where in 
the subscript of the susceptibility, C refers to the vari- 
able conjugated to the field T: C — j^, H being the 
Hamiltonian of the perturbed system. The FDT states 
that the susceptibility XAc{t) is related to the equilib- 
rium correlation function CAc(t) — (A(t)C(0)) through: 
J XAc( s )ds = -j^fCAc{t) where T is the thermal bath 
temperature, and the correlation function is estimated 
at equilibrium. A sensitive way of probing the devia- 
tion from this relation in nonequilibrium systems, which 
has been extensively used for example in glassy systems 
[15, 16] consists in determining separately the integrated 
susceptibility function and the correlation function, and 
in plotting them in a parametric plot with the time as pa- 
rameter. The slope of the curve is then interpreted as the 
inverse of an effective temperature, which may depend on 
the time scale [15]. 

For the system under study, we obtain the integrated 
response to an external force F by applying the force in a 
stationary configuration at t = 0, and following the evo- 
lution of the particle center of mass velocity U(t). The 
parametric plot involves then the average velocity divided 
by the applied force, fi(t) = (U(t))/F = L xux(s)ds, 
versus the integrated velocity auto correlation function 
Cux{t) — J Q Cuu(s)ds = T>(t). To obtain the response 
function from the ensemble averaged particle velocity 
(£/(£)), we have run simulations starting from 1000 in- 
dependent configurations of the system and tracked the 




FIG. 2. (a) Integrated velocity autocorrelation func- 
tions of the particle (from bottom to top: ksTp/e = 
0.75,1.5,2,2.5,3,3.5). (b) Einstein temperature Te and ki- 
netic temperature Tk as a function of the particle tempera- 
ture T p ; the contact temperature T c is also plotted for com- 
parison. Lines are guides for the eye. When not indicated, 
uncertainties are below the symbol size. 
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FIG. 3. Integrated response function as a function of the in- 
tegrated VACF of the nanoparticle, for fcaTp/e = 0.75 (equi- 
librium) and 3.5. Inset- Temperatures extracted from the fit 
of the main graph's curves at small and large times, as a func- 
tion of the particle temperature, note that the lines are not 
merely guides to the eye, but correspond to the data deter- 
mined independently and already reported in Fig. 2 for the 
kinetic and Einstein temperature. 



position of the Brownian particle before a steady state is 
attained (corresponding to times smaller than t c ). This 
enabled us to obtain good statistics for the ensemble av- 
eraged velocity of the particle, in particular during the 
early stage of the transient t <C i c - 

Figure 3 shows the resulting response/correlation para- 
metric plot, for the different temperatures considered. 
When T p = To, the nanoparticle is at equilibrium before 
the external force is applied, and the fluctuation dissi- 
pation theorem is obeyed. For higher values of the par- 
ticle temperature T p , the velocity (U(t)) depends non 
linearly on the integrated VACF and the fluctuation dis- 
sipation ratio is time dependent. This is particularly vis- 
ible for the highest temperature considered in Fig. 3 



3.5e/£;b, where the two slopes 4^ at small and 



large V differ markedly. From these two slopes, it is pos- 
sible to define two temperatures T^ and T E characteriz- 
ing the response of the system respectively at short times 
and long times. The inset of Fig. 3 compares these two 
temperatures to the kinetic and Einstein temperatures 
defined before. Strikingly the short time effective tem- 
perature T K is very close to the kinetic temperature of 
the nanoparticle Tk, while the long time effective temper- 
ature T E is close to the Einstein temperature Te. There- 
fore our system, in spite of its simplicity, exhibits a "two 
temperatures" behavior on the two different time scales 
that are separated by the typical scale set by the loss of 
memory in the initial velocity. The short time, fast tem- 
perature sets the kinetic energy of the particles, while 
the Einstein temperature which probes the steady state 
response is determined by the long time behavior of the 
integrated response. 

For a system in contact with a thermal bath and driven 
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FIG. 4. (a) Transient fluctuation temperature Tt = 
{SWt) /2(Wt) as a function of the time t, for different tem- 
peratures T p of the nanoparticle. From bottom to top: 
k B T p /e = 0.75,1,1.5,2,2.5,3,3.5. (b) Transient fluctuation 
temperature Ttft obtained with the long time limit of T t as 
a function of the particle temperature T p . The lines corre- 
spond to the data for the kinetic and Einstein temperature in 
Fig. 2. 



245 out of equilibrium, the bath temperature plays also a key 

246 role in quantifying the fluctuations of the work from an 

247 external forcing [10]. Two situations have to be distin- 

248 guished depending on the time window analyzed. If we 

249 follow the evolution of a system in the transient regime 

250 before a steady state is reached, starting from a system 

251 at equilibrium, the transient fluctuation theorem (TFT) 

252 predicts: 



P(W t )/P(-Wt) = exp(W t /fc B T), 



(3) 



where P(Wt) is the density probability of the work Wt- 
In this equation Wt is the work from the external force 
F, i.e. W t = J Q U(s)Fds and T is the temperature 
of the thermal bath. On the other hand, in a a sta- 
tionary situation, the steady state fluctuation theorem 
(SSFT) predicts P{W t )/P{-W t ) exp(W t /fc B T) when 
t 3> t c where t c denotes a typical equilibrium correla- 
tion time. In the SSFT, the work Wt is estimated along 
a trajectory of length t: W t — f t %+ U(s)Fds, where 
an average on different values of the initial U may be 
performed. We have tested these fluctuation relations 
for the heated Brownian particles, again applying an 
external force F = lOe/er at t — and recording the 
statistics of the work using 1000 independent configura- 
tions. It turned out however that the distribution of the 
work Wt was too noisy to determine accurately the ratio 
P(Wt)/P(—Wt) and critically assess the validity of the 
fluctuation theorems discussed above. To extract an ef- 
fective temperature measuring the fluctuations of W t , we 
have used the observation that the statistics of the work 
Wt is to a good approximation Gaussian. Under these 
conditions, it is trivial to show that the distribution of 
Wt obeys a law similar to Eq. 3 with an effective tem- 
perature T t = (5Wt)/2(Wt). Note that strictly speaking 



the TFT implies that T t = T is independent of t. In Fig. 
4. a we have shown the evolution of T t as a function of 
the time t for different temperatures T p of the nanopar- 
ticle. For all the temperatures considered, the initially 
small values of (W t ) leads to a large uncertainty in the 
value of T t . For longer times t > 5r, the temperature 
T t is approximately independent of the time t. We will 
denote Ttft(T p ) the value of the effective temperature 
Tt in this regime. Figure 4.b displays the evolution of 
Ttft as a function of the temperature of the nanoparti- 
cle T p . It is clear that the resulting Ttft is very close 
to the kinetic temperature Tk characterizing the particle 
dynamics on short time scales. While we are not aware 
of a theoretical analysis of this situation, we believe the 
reason for this proximity lies in the fact that the main 
contribution to fluctuations in the work function corre- 
sponds to the time regime in which the velocity is still 
correlated to its value at t = 0, i.e. the same time regime 
in which the fluctuation dissipation ratio corresponds to 
the "fast" temperature. 

Our work shows that, even in a conceptually rather 
simple system, in a nonequilibrium steady state, a de- 
scription in terms of a Langevin model involving a single 
temperature is far from trivial. Further generalization 
and interpretation of the behavior of interacting parti- 
cles in terms of Langevin models and a single noise tem- 
perature is expected to suffer similar difficulties, as can 
already be inferred from the results of [2]. It would be 
interesting to explore, if the recent extensions of fluctu- 
ation dissipation theorems proposed in refs [17, 18] can 
be applied to the present case, i.e. to identify observ- 
ables for which a response-correlation proportionality re- 
lation holds. Even so, the resulting observables are likely 
to be different from those that are naturally measured 
in experiments or simulations. We also note that, with 
the present observables, experiments using optical tweez- 
ers with a strongly absorbing particle could be used to 
probe the different temperatures investigated here, with 
the exception of the kinetic one. We expect that such 
experiments will be able to detect a deviation from equi- 
librium of the order of magnitude reported here. 
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